function meanerr(time,a,m,~,sig0,sig)

load('orthbase');

%exactmean=zeros(1,100);

exactmean=sig0*exp(-a*time)+m*(1-exp(-a*time));
hold on
plot(time,exactmean);
plot(time,sig(:,1),'r-');
title('Mean')
legend('exact mean','wce mean')
